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ABSTRACT 

We present updated estimates of Galactic foreground emission using seven 
years of WMAP data. Using the power spectrum of differences between multi- 
frequency template-cleaned maps, we find no evidence for foreground contami- 
nation outside of the updated (KQ85y7) foreground mask. We place a 15 fiK 
upper bound on rms foreground contamination in the cleaned maps used for cos- 
mological analysis. Further, the cleaning process requires only three power-law 
foregrounds outside of the mask. We find no evidence for polarized foregrounds 
beyond those from soft (steep-spectrum) synchrotron and thermal dust emission; 
in particular we find no indication in the polarization data of an extra "haze" 
of hard synchrotron emission from energetic electrons near the Galactic center. 
We provide an updated map of the cosmic microwave background (CMB) using 
the internal linear combination method, updated foreground masks, and updates 
to point source catalogs using two different techniques. With additional years 
of data, we now detect 471 point sources using a five-band technique and 417 
sources using a three-band CMB-free technique. In total there are 62 newly de- 
tected point sources, a 12% increase over the five-year release. Also new are tests 
of the Markov chain Monte Carlo foreground fitting procedure against systemat- 
ics in the time-stream data, and tests against the observed beam asymmetry. 

Within a few degrees of the Galactic plane, the behavior in total intensity 
of low-frequency foregrounds is complicated and not completely understood. 
WMAP data show a rapidly steepening spectrum from 20-40 GHz, which may be 
due to emission from spinning dust grains, steepening synchrotron, or other ef- 
fects. Comparisons are made to a 1-degree 408 MHz map (Haslam et al.) and the 
11-degree ARCADE 2 data (Singal et al.). We find that spinning dust or steepen- 
ing synchrotron models fit the combination of WMAP and 408 MHz data equally 
well. ARCADE data appear inconsistent with the steepening synchrotron model, 
and consistent with the spinning dust model, though some discrepancies remain 
regarding the relative strength of spinning dust emission. More high-resolution 
data in the 10-40 GHz range would shed much fight on these issues. 

Subject headings: cosmic background radiation — cosmology: observations — 
diffuse radiation — Galaxy: halo — Galaxy: structure — ISM: structure 
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Introduction 



The Wilkinson Microwave Anisotropy Probe (WMAP) was launched in 2001 to observe 
the cosmic microwave background (CMB). In addition to measuring the CMB, WMAP, 
like any full-sky CMB experiment, also observes emission from our own Galaxy. With five 
frequency bands centered at 23, 33, 41, 61, and 94 GHz (respectively denoted as K, Ka, 
Q, V, and W bands), full sky coverage, polarization sensitivity, and control of systematics 
to the sub-percent level, WMAP is able to measure diffuse (1° and larger) emission with 
precise temperature calibration. In this paper we analyze seven years of WMAP data in 
order to better characterize Galactic foreground emission, the removal of which will be one 



of the largest challenges to future CMB experiments (Dunkley et al. 2009b) 



This paper is part of a suite of papers describing the full details of the WMAP seven-year 



data release. An overall description of sky maps and basic results is in Jarosik et al. (2010) 



which also includes a description of the beam modeling used to produce maps smoothed to 
the common resolution of a 1° FWHM Gaussian. These maps serve as the starting point 



for foreground analysis in this work. Larson et al. (2010) describe the generation of power 



spectra from CMB maps, and Komatsu et al. (2010) discuss the cosmological implications of 



the spectra. Weiland et al. (2010) detail measurements of celestial calibrators, and Bennett 



et al. (2010) investigate the status of some potential anomalies found in WMAP data. 



The layout of this paper is as follows. Updates to masks and foreground fitting processes 
are described in Section [2] A comparison of WMAP data to that recently taken by the 



ARCADE instrument ( Singal et al. 2009 ) is discussed in Section 3 Results of the fits and 



their implications for specific foreground emission processes are discussed in Section |4} A 
discussion of systematics follows in Section |5| Point sources and an update to the point 
source catalog are found in Section [6j Lastly, conclusions can be found in Section [7| 



1.1. Science Overview 



There are three primary mechanisms for diffuse Galactic radio emission. Relativistic 
electrons interact with the Galactic magnetic field to produce synchrotron emission, for which 



the standard template is 408 MHz data compiled by Haslam et al. (1981). Less energetic 



electrons scatter from each other and ionized nuclei to produce free-free radiation (also 



known as thermal Bremsstrahlung) , which can be traced with Ha line emission (Finkbeiner 



2003). Finally, dust grains emit a modified blackbody spectrum through excitation of their 



vibrational modes, for which the standard template is the fit of Finkbeiner et al. (1999) to 



data from the Infrared Astronomical Satellite (IRAS) and the Cosmic Background Explorer 
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( COBE). Dust grains may also emit radiation through rotational modes or other excitations 
dPraine fc Lazarian|[l998a||b[ p99| . 



WMAP was designed to measure near the frequency where the ratio of the CMB 
anisotropy to the rms fluctuations of all three foregrounds is at its maximum, to minimize 
foreground contamination. This also implies that two or more foreground components will 
be of comparable amplitude and that they will be relatively weak. Foreground templates, 
however, are best made by observing a foreground process at a frequency where it dominates 
the total emission. Hence there will always be some extrapolation involved when attempting 
to account for foregrounds on top of CMB observations. 

So how well does the extrapolation work? Simple power-law extrapolation of the 408 



MHz synchrotron template from Haslam et al. (1981) does not explain very much of the 



observed emission at 20-40 GHz. Whether this is due to a new low-frequency emission 
process, errors in the extrapolation due to spatial variation in the spectral index, or both, is 



difficult to determine. Targeted observations of individual regions (Scaife et al. 2009; Tibbs 



et al.||20l"0 ) suggest a spinning dust-like component, but a model consistent across size scales 
and data sets remains elusive. 

Free-free emission is extrapolated from a map of Ha, corrected for dust extinction using 
a reddening map based on 100 yum data ( Schlegel et al.||1998 Bennett et al.|200^ ). Variations 
in electron temperature cause some uncertainty in this extrapolation, but the larger effect 
is likely uncertainty in the reddening correction. The overall ratio of radio to Ha brightness 



comes out lower than expected (Bennett et al. 2003); nevertheless, the template otherwise 



matches quite well with observations at 30-60 GHz, where free-free emission from the Galactic 
disk is particularly dominant. 

The dust extrapolation has so far been tested least precisely by CMB experiments. 



While the model of Finkbeiner et al. (1999) incorporates COBE FIRAS data all the way 



down to 60 GHz, the uncertainty at those frequencies is large; most of the dust model 
comes from information at 100 and 240 /xm. While the spectral index of dust at frequencies 
below 300 GHz has not yet been measured to enough accuracy to challenge the model, the 
morphology matches observations at lower frequencies, though some experiments suggest 



overall brightness levels different from the predictions (Veneziani et al. 2010; Culverhouse 



et al. 2010). 



Analysis of data from previous WMAP releases has shown that CMB maps from different 



foreground removal techniques agree to within 11 /xK (Gold et al. 2009) on average in the 



low Galactic emission regions used for CMB anisotropy measurements, though this does not 
provide an absolute limit to the amount of contamination. Even when templates are not 
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directly used for foreground removal, they provide an important guide for the construction 
of masks and other foreground cleaning methods. 

Several systems of units are in use throughout this work. Point sources are reported 
in flux units (Jansky), where the power-law index is denoted by a such that flux follows 
S ~ z/°. Foreground modeling is most easily done in units of antenna temperature, defined 
by using the Rayleigh- Jeans limit of a blackbody spectrum (for which S ~ z/^) to convert 
flux per solid angle to a temperature. In these units the power-law index is denoted as 
/3 = a — 2. y\rM.AP'''s, frequency range is not quite in the Rayleigh- Jeans limit for a 2.7 
K blackbody, so there is a frequency-dependent conversion factor a(z/) = (e^ — l)^/x^e^ 
(where x = hv /kBTcmh) to convert antenna temperatures to thermodynamic temperatures 
convenient for CMB analysis. 



2. Seven-year Foreground Fits 
2.1. Masks 

Foreground removal always has some uncertainty, so it is useful to mask part of the sky 
where foregrounds are too bright for CMB analysis. As in the five-year analysis, the starting 
points for the masks are K- and Q-band-average maps smoothed to one-degree resolution. 
The maps are then converted to foreground-only maps by subtracting off an estimate of the 



CMB using the internal linear combination (ILC) method (see Hinshaw et al. 2007, and 



Section 2.2). A cumulative histogram in each band is formed to find the flux level above 
which a given percentage of sky can be cut, and the union of the pixels cut from each band 
at a given flux level is used to define a mask. We used two masks for most further analysis, 
based on cuts which leave 75% and 85% of the sky; these are denoted as KQ75 and KQ85, 
respectively. 

For the seven-year analysis, the diffuse foreground masks have been extended based on 
a analysis of residuals after foreground subtraction. Starting with foreground-reduced 
maps, differences are taken between bands {Q — V and V — W m thermodynamic units), 
eliminating any CMB signal. Ideally the only thing left in the resulting maps would be noise; 
in practice there are visible residuals near the Galactic plane. Given the noise per pixel of 
the maps, it is possible to compute a map of the for each pixel. 



After degradation to HEALPix N^i^^ = 32 (see Gorski et al. 2005 for a description of 
this pixelization scheme), regions of 4 or more contiguous pixels with higher than four 
times that of the polar caps are identified and used to define two new masks, one from each 



difference map. These are then combined with the previous KQ75 and KQ85 masks (Gold 
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et al. 2009) used for the five-year analysis. After promotion back to full resolution, edges of 
the mask are smoothed with a 3° FWHM Gaussian. The resulting changes to the final mask 
are primarily around the edge of the Galactic cut, particularly in the Gum and Ophiuchus 
regions. The additional sky fraction cut from the KQ85 masked sky is 3.4%, and from the 
KQ75 masked sky is 1.0%. 

These expanded masks are then combined with the point source mask as in previous 
releases, which has been updated with newly detected sources. Also, point sources brighter 
than 5 Jy have had the radius of their cut extended from 0°.6 to 1°.2, in order to minimize 
confusion at low frequencies where the instrument beam is large. The new masks, which 
we denote KQ75y7 and KQ85y7, are shown in Figure [TJ and are available on the LAMBDA 
Web sit(j5 In total 70.6% (KQ75y7) and 78.3% (KQ85y7) of the sky now remains after the 



masking process. 



2.2. Internal Linear Combination Method 



The ILC method implemented by WMAP is a technique largely blind to assumptions 
about the frequency spectrum of foreground emission, which produces CMB maps with little 
visible foreground contamination. The ILC is a weighted combination formed from all five 
frequency bands, which are smoothed to a common 1° FWHM Gaussian beam using the 



symmetrized beam window functions produced by the beam analysis (Jarosik et al. 2010). 



The coefficients used to weight each individual frequency band are those that minimize the 
variance of the resulting map under the constraint that the sum of the coefficients is unity, 
which ensures that the CMB portion of the signal is passed through unaltered. 

The details of the algorithm used to compute the WMAP seven-year ILC map are the 
same as that described in the three- year analysis (Hinshaw et al. 2007). In particular, we 



perform a bias correction step which uses simulations to estimate and correct for the tendency 
of the ILC method to produce CMB maps anti-correlated with foreground fluctuations (for an 
overview of potential ILC pitfalls see Vio fc Andreani|2009 ). We have found this technique to 
be robust when applied to WMAP data: the variance between the ILC map and CMB maps 



made with other techniques is less than 116 /iK (Gold et al. 2009). Similar techniques by 



other authors have given CMB maps consistent with WMAP^s best-flt cosmological results 



(Kim et al. 2008). 



Rather than use a single set of coefficients for the whole sky, to allow for variations in 



^http : / /Icunbda . gsf c . nasa . gov/ 
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Galactic composition we subdivide the sky into 12 regions and find the ILC coefficients for 
each, shown in Table [TJ All but region lie along the Galactic plane. We retain the same 
number of regional subdivisions of the sky and their spatial boundaries remain unchanged 

. The frequency weights for each 
region are slightly different, however, reflecting the most recent updates to the calibration 
and beams. Figure |2] shows the difference between the seven-year and five-year ILC maps, 
which is dominated by a small change in the dipole. The seven-year ILC map is available on 
the LAMBDA Web site. 



from the previous years (for details see [Hinshaw et aL][2007| 



2.3. Maximum Entropy Method 



The maximum entropy method (MEM) is a spatial and spectral fit that uses external 
templates, intended to distinguish between different emission sources. By design, the MEM 
output tends to revert to these templates in regions of low signal-to-noise. Thus the MEM 
results are most interesting in regions with higher signal. 



The seven-year MEM analysis is largely unchanged from previous work (Hinshaw et al. 



2007 Gold et al. 2009). As before, the analysis is done in all bands on the 1° smoothed sky 



maps, with the ILC map subtracted. The zero level of each map is set such that a esc \ b\ 
fit, for HEALPix A^side = 512 pixels at 6 < —15° and outside of the KQ85y7 mask, yields a 
value of zero for the intercept. The maps are degraded to HEALPix A^sidc = 128 pixelization, 
and a model is fit for each pixel p. Rather than simply minimize x^, the MEM minimizes a 
function 



^ = x' + A5^T,(p)ln 

c 



Tcjp) 
[ePc{p)\ 



T. 



Here Tc and Pc are the model brightness and template prior for foreground component c 
(e is the base of natural logarithms). The second term is what enforces the prior when 
the signal-to-noise becomes low, and the parameter A sets the threshold for the transition 
from signal-dominated to noise-dominated behavior. The spectra of the free-free and dust 
components are fixed power laws, with /3 = —2.14 for free-free and /3 = +2.0 for dust. An 
iterative procedure uses residuals from the fit at each iteration to adjust the spectrum of 
the synchrotron component for each pixel. Hence any anomalous component such as electric 
dipole emission from spinning dust is included in the synchrotron component. The adopted 



priors are unchanged from previous analyses and are based on the 408 MHz map of Haslam 



et al. (1981) with an extragalactic brightness of 5.9 K subtracted (Lawson et al. 1987) for 



synchrotron, an extinction-corrected Ha map (Finkbeiner 2003 Bennett et al. 2003) for 



free-free, and model 8 of Finkbeiner et al. ( 1999 ) for dust 
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The prior map and output map are shown in Figure |3] for each foreground component. 
The zero level of the output synchrotron map is slightly lower (~ 50 fiK) than that of the 
synchrotron prior. This reflects the difference between the zero level of the K-hand esc |6| 
normalized map and that of the prior. For comparison, the la uncertainty in the prior zero 
level, based on the quoted uncertainty in the 408 MHz map zero level, is 27 fiK. Also, there 
is a dependence on the adopted extragalactic brightness at 408 MHz. If the ARCADE 2 



value (Fixsen et al. 2009) were used, the zero level of the prior would be ~ 37 jjK below 



that of the esc \ b\ normalized map. Figure [s] can be compared with Figure 5 of (Hinshaw 



et al. 2007) to see the improvement in signal-to-noise ratio of the output maps between the 



three-year and seven-year analyses. 

Differences between seven-year MEM maps and five-year MEM maps are shown in 
Figure |4j The seven-year MEM foreground component maps tend to be slightly brighter 
than the five-year versions at mid to high northern Galactic latitudes. This is due to small 
dipole differences between the seven-year and five-year sky maps, which are caused by a 
combination of a change in the calibration dipole and small (less than 0.2%) changes in 
radiometer calibrations between the seven-year and five-year analyses. The seven-year and 
five-year foreground component maps are in better agreement at southern Galactic latitudes 
because this is where zero level normalization of the sky maps is determined by esc \ b\ fitting. 
The MEM maps are available on the LAMBDA Web site. 



2.4. Template Cleaning 

WMAP continues to use a template cleaning method to produce the foreground-reduced 



maps used for power spectrum analysis (Hinshaw et al. 2007 Page et al. 2007). For tem- 



perature maps, the templates are a K— Ka difference map, an extinction-corrected Ha map. 



and a dust map Finkbeiner et al. (1999). For polarization, the templates are the K-hand 



map for synchrotron, and a dust model described in detail below. 

The temperature cleaning is applied to seven-year Q-, V-, and VT-band maps {K and 
Ka are used for a template). The model has the form 



(2) 



where p indicates the pixel, the frequency dependence is entirely contained in the coefficients 
bi, and the spatial templates are the WMAP K-Ka temperature difference map [Tk — Txa), 
the Finkbeiner (2003) composite Ha map with an extinction correction applied (/hq), and the 



Finkbeiner et al. (1999) dust model evaluated at 94 GHz (Mjust)- Because the first template 



has contributions from both synchrotron and free-free emission, foreground parameters are 
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a mixture of and &2('^)- For free-free emission, the ratio of K-hand radio temperature 
to Ha intensity is 

5ff(z/) - 0.552 6i(z/) ^ ' 

where 5'ff(z/) is the free- free emission spectrum converted to thermodynamic temperature 
units, normahzed to unity at i^-band, and is assumed to be a power-law in antenna tem- 
perature with /3 = —2.14. The synchrotron spectral index (relative to i^-band) is found 
via 

^ \og[0.67b,{u)a{u)] 

where a(z/) is the conversion factor from antenna temperature to thermodynamic units. 

The coefficients of the model fit to the seven-year data are presented in Table [2] Small 
changes in the seven-year coefficients compared to previous values reflect small changes in 
the absolute calibration and beam proflles. 

For polarization cleaning the maps are degraded to low resolution {Ngi^e = 16). The 
model has the form 

[Q(z/,p), f/(z/,p)] model = ai{u)[Q{p),U{p)]K + a2{u)[Q{p),U{p)]dust (5) 

The templates used are the WMAP i^-band polarization for synchrotron {[Q,U]k), and a 
low resolution version of the dust template used above with polarization direction derived 
from starlight measurements ([Q, t^jdust) and a geometric suppression factor to account for 



the magnetic fleld geometry (Page et al. 2007). The coefficients of the model fit to the 



seven-year data are in Table |3j For polarization, the template maps are assumed to have a 
one-to-one correspondence with foreground emission, so the spectral indices for synchrotron 
and dust are simply the power-law slopes of the coefficients ai(z/) and a2{i'). If the dust 
model is correct then the ratio a2/&3 gives the polarization fraction; for W^-band this is 
-6%. 

The full- resolution {Nside = 512) foreground-reduced Stokes Q and U maps were pro- 
duced using the same cleaning coefficients as derived for the low-resolution maps, but with 
full-resolution templates. The i^-band and dust intensity templates can be produced at 
full resolution from available data, and the starlight polarization map used to determine 
polarization direction was upgraded to full resolution using nearest-neighbor sampling. The 
templates subtracted from WMAP data are smoothed to 1° FWHM, potentially leaving 
artifacts in the foreground-reduced maps due to small-scale power or beam asymmetries. In 
practice, we find no sign of these effects, as discussed in Section 4.1 and Section [5j All data 
sets used for templates are available on the LAMBDA Web site. 
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2.5. Markov Chain Monte Carlo Fitting 

We again perform a pixel-based Markov chain Monte Carlo (MCMC) fitting technique to 



the five bands of WMAP data. Our method is similar to that of Eriksen et al. (2007), but we 
focus more on Galactic foregrounds rather than CMB. The fit results of the five-year release 
have been reproduced, with the "base" model, which uses three power-law foregrounds, 
producing virtually the same reduced P^^^ pixel. The MCMC fitting has benefited from 
further understanding of the zero point of the maps. We have used the 408 MHz map of 



Haslam et al. (1981) with a zero-point determined using the same esc |6| method as for the 



WMAP data, and investigated the effect on the fit of error in this zero-point. 

The MCMC fit is performed on one-degree smoothed maps downgraded to HEALPix 
-^side = 64. A MCMC chain is run for each pixel, where the basic model is 

/ u \ ^'^'''^ ( v\^^ ( vV" 

^'"'"^•U) +aMT.„. + r.(-) (6) 

for the antenna temperature. The subscripts s, /, d stand for synchrotron, free-free, and dust 
emission, vk and vy^ are the effective frequencies for K- and ly-bands (22.5 and 93.5 GHz), 
and aiv) accounts for the slight frequency dependence of a 2.725 K blackbody using the 



thermodynamic to antenna temperature conversion factors found in Bennett et al. (2003). 
The fit always includes polarization data as well, where the model is 

Q{^) = QA — \ +Qd — +a(i.)Qcmb (7) 



U[v) 



Vk ) \ J 



Us( — ] +UA + a Z/ t/cmb 



for Stokes Q and U parameters. Thus there are a total of 15 pieces of data for each pixel 
(T, Q, and U for five bands). 

As for the five-year release, the noise for each pixel at A^side = 64 is computed from maps 
of Aobs at Aside = 512. To account for the smoothing process, the noise is then rescaled by 
a factor calculated from simulated noise maps for each frequency band. The MCMC fit 
treats pixels as independent, and does not use pixel-pixel covariance, which leads to small 
correlations in between neighboring pixels. This has negligible effect on results as long as 
goodness-of-fit is averaged over large enough regions. 

We fit three categories of models. All use A'-band as a template for the polarization angle 
of synchrotron and dust emission, so Us and Ud are not independent parameters, identical 
to the previous analysis. All models also fix the free-free spectral index to /3/ = —2.16, a 
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slight change from /3/ = —2.14 used in the previous analysis. This change was motivated as 
an attempt to better match the effective spectral index at Q- and l^-bands, due to their use 
in cosmological analysis, but was not found to make a difference to the fits. 

The "base" model uses three power-law foregrounds, where the synchrotron spectral 
index Psiy) is taken to be independent of frequency but may vary spatially, and the dust 
spectral index Pci is allowed to vary spatially. We assume the same spectral indices for 
polarized synchrotron and dust emission as for total intensity emission. This model has a 
total of 10 free parameters per pixel: T^, T/, T^, Tcmb, A, /3d, Qs, Qd, Qcmb, and Ucmh- 

A steepening synchrotron model uses the same three foregrounds but allows for a steep- 
ening of the synchrotron spectral index by adding a new parameter c^, defined by 



/5^W= . , „ ..... ■ (9) 




For the steepening model the dust spectral index is fixec^^ to /3a = +2.0. Therefore this 
model also has 10 free parameters per pixel. 

For models with a spinning dust component, another term is added to equation [6] 

Tsdii^) = ^^1^^^" (10) 
exp(z//i/sd) - 1 

The spinning dust component is assumed to have negligible polarization, as theoretical ex- 



pectations for the polarization fraction are low compared to synchrotron radiation (Lazarian 



fc Draine||2000 ), and the polarization data thus far show no evidence that such a component 



is necessary (see Section 4.5). The spinning dust amplitude Ag, was allowed to vary spatially 
as a new parameter. Both and [3, were fixed to —3.0 and +2.0, respectively, to avoid 
degeneracies from having too many parameters in the fit. Allowing Vgd to spatially vary was 
not found to result in any improvement of the fit, but fits were performed with different 
global values of Vgd to find the best overall value. Thus with fixing of the spectral indices, 
this model has nine free parameters per pixel. 

MCMC fits for the seven-year release were performed with the addition of the 408 MHz 



data compiled by Haslam et al. (1981). The error on the zero point for this data was es- 



timated in that work to be ±3 K, with an overall calibration error of 10%. Lawson et al. 



^^The precise choice of dust index here and for the spinning dust model does not make much difference; 
when allowed to vary it is poorly constrained by the MCMC fits and uncorrelated with the synchrotron or 
free- free components ( Gold et al.||2009 ). 
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(1987) use a comparison with 404 MHz data to find a uniform (presumably extragalactic) 
component witfi a brigfitness of 5.9 K. As tfie MCMC metliod treats all input maps equally, 
for consistency we estimate and subtract off a nominal zero point offset of 7.4 K, as de- 
termined by the same esc |6| method we use for the V\[MAF sky maps. However, the 408 
MHz data resembles a esc |6| behavior much less than the MIMAP data, due to the increased 
relative prominence of large-scale features such as the Northern Polar Spur. Therefore we 
attempted the esc |6| fitting procedure on different hemispheres and with different cuts, and 
estimate the uncertainty in procedure to be ±4 K. MCMC fits were run for each model with 
zero points of 3.4 K and 11.4 K in addition to the nominal value, and the effect of these on 
foregrounds is discussed in Section 4^ A full set of maps and MCMC variance estimates 
for the three models is available on the LAMBDA Web site. 



3. Comparison with ARCADE 2 

The ARCADE collaboration has made available absolute temperature measurements of 



Galactic emission for part of the sky (Kogut et al. 2009). ARCADE observations do not 



cover the full sky and the instrument's beam is significantly larger than WMAP''^. Therefore 
we limit our comparison to two regions where ARCADE's scan crosses the Galactic plane 
and observes the brightest emission, the first at Galactic longitude of 34° and the second at 
93°. 

Figure [5] shows the Galactic spectrum for these two regions. WMAP data have been 
smoothed to 11.6°, to match the ARCADE resolution. The ARCADE maps have had the 
CMB monopole removed, and the WMAP maps have had CMB anisotropics removed using 
the ILC map (though this has little effect). The ARCADE data have not had any extra- 



galactic component (as found by Fixsen et al. 2009) removed. Instead, all maps have been 



treated as equally as possible, removing a zero-point by fitting a esc \b\ model to the available 
data and subtracting the constant term. 

The uncertainty in this zero-point subtraction is largest for ARCADE due to the limited 
sky coverage of the experiment. We tested the zero-point subtraction by fitting to several 
partial-sky subsets of the full-sky V\fMAP maps, and find that the variations imply an 
uncertainty in the ARCADE points of up to 15% of the CMB-subtracted flux. Also included 
is the 408 MHz map as a reference point at low frequency. As discussed in the previous 
section, the esc |6| model performs most poorly for this map, with uncertainties of ±4 K. 
However, in these two regions the emission is bright enough that this is still less than 10% 
of the total emission. 
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Two fits were applied to the data in each region. The first used three power-law fore- 
grounds: synchrotron, free-free, and dust, where the spectral indices for synchrotron and 
dust were left free for the fit. The second fit added a spinning dust component using the 
functional form of Equation 10, with the amplitude and Ugd as free parameters. Because 
the maps are highly smoothed, errors are dominated by systematic issues and difficult to 
characterize. We chose to use 2% fractional error for WMAP and 5% fractional error for 
other observations when performing the fit. Using larger errors does not remove the sharp 
difference in between the two models unless the errors are taken to be larger than 50%. 
For the fit without spinning dust, the ARCADE data were not used in the fit as they were 
found to be incompatible with such a model. 

The resulting fits are shown in Figure [5} with the spinning dust fit in blue and the 
power-law-only fit in red. The top panels show the data and fits in absolute temperature 
units after monopole subtraction. The bottom panels show the same data and fits, but where 
all temperatures have had the 0.408-22 GHz slope divided out, to facilitate comparison with 
Figure 9 of |Kogut et~al] ( |2009D . The ARCADE data show a clear deficit over the 3-10 
GHz range, which cannot be explained with power-law foregrounds alone; a fit including a 
spinning dust component is much more consistent. Dotted, dashed, and dash-dotted lines 
in the figure show the contribution of each individual component to the total, with thermal 
dust and spinning dust shown together. In the spinning dust model, synchrotron emission 
is weak in the WMAP bands, where free-free is the dominant emission process. At 93° 
longitude the spinning dust emission is approximately as bright as the free-free emission at 
23 GHz, and at 34° longitude it is several times fainter at all frequencies. 



4. Foreground Results 

4.1. Residuals in Template-Cleaned Maps 

As a test of the template-based foreground subtraction process, power spectra of dif- 
ference maps were made. Figure |6] shows the power spectrum of the difference between the 
foreground-reduced Q-band and W^-band maps, with the point source contribution to the 
power spectrum subtracted off. Averaging over bins of A£ = 50, no bin with more than 120 
fiK^ of power is seen, with an upper hmit of ~ 220 /iK^ in power (15 /xK in amplitude), and 
the results are consistent with zero within the expected error. For comparison, CMB power 



in the range 30 < £ < 500 is 1000 /iK^ or more (Larson et al. 2010). Differences between 
foreground-reduced V^-band and VT-band were also computed, and the power in that case 
was even smaller. 
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4.2. Polarization Power Spectra of Synchrotron and Dust 

While Galactic foregrounds are not fully described by a two-point function (i.e. an 
angular power spectrum), due to the importance of the CMB it is often useful to examine 
the power spectrum of foregrounds. Specifically, the relevant quantity to calculate is the 
contribution of foreground emission to the angular power spectrum in a particular patch of 
sky of interest for CMB analysis. 

A general trend of ^{^ + 1)C^ ~ was found from examination of raw polarization 
data outside the P06 mask ( Page et aL]|2007 ), as the result of a combined fit to WMAP data 



in both multipole and frequency space. With the MCMC fitting procedure it is possible 
to separate polarized synchrotron from dust and examine the two components individually. 
The results, shown in Figure [7], show behavior largely consistent with the previous analysis. 

In detail, MCMC maps from the "base" model including Haslam data were used. Power 
spectra from the spinning dust model MCMC maps were also inspected and found to be 
nearly identical at large scales. A union of the polarization analysis mask and the mask of 
pixels fiagged by the MCMC was applied, and the Cf^ and Cf^ spectra were computed 
for both synchrotron and dust. As the MCMC process uses one-degree smoothed maps, an 
appropriate correction for the beam window function was applied. Each power spectrum 
was then fit with a model consisting of a power-law plus a pixel noise term 

£(£ + l)Cf ^/27r = BJ"^ + i{i + l)N^, (11) 

where Be is the amplitude for foreground component c, m is the power-law index, and N the 
noise amplitude. 

Values for the fit parameters and an estimate of their errors can be found in Table |4j 
Because the power spectra are taken from highly processed maps, detailed error propagation 
is difficult. We used the diagonal portion of the published Ci Fisher errors plus cosmic 
variance to perform the fit; covariance between multipoles will cause the true errors to be 
somewhat larger. If all foreground power spectra are assumed to have the same power-law 
behavior, then the weighted mean m = —0.67 ± 0.24. 



4.3. Free-free Emission 

That the ratio of radio brightness to Ha intensity from the MEM fits is consistently 
lower than the expected value has long been of concern. The MCMC fits offer some insight, 
though unfortunately do not resolve the difference. The most important difference between 
the MEM and the MCMC fits in this case is that the MEM uses the Ha template as a prior 
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in low signal-to-noise regions, while the MCMC fit does not. The result is that in regions of 
low signal, the degeneracy between synchrotron and free-free causes the MCMC uncertainty 
in free-free brightness to be large enough to accommodate a large range of possible radio 
to Ha ratios. Therefore it becomes necessary to exclude low signal-to-noise regions when 
calculating the ratio from the MCMC maps. 

Due to the uncertainty in the reddening correction to the Ha map itself, it is also 
customary to exclude regions where the Ha optical depth due to reddening is greater than 
some value. This unfortunately excludes regions that would otherwise have high signal-to- 
noise. These two cuts together exclude much of the MCMC maps as unsuitable for analysis. 
The remaining portion of sky contains bright, mostly discontiguous free- free regions which 
are also low on dust (and therefore Ha extinction). The largest of these is a region around 
Gum nebula. 

Starting with the free-free maps made from the MCMC process, we define a signal-to- 
noise ratio (SNR) map as the free-free amplitude divided by the square root of the MCMC 
variance. We then keep only pixels with SNR > 10, r < 1, and no MCMC error fiags. The 
pixels that remain are largely concentrated in three regions, the Gum nebula, the Ophiuchus 
complex, and the Orion/Eridanus bubble. The Gum region contains nearly half of the pixels 
surviving the cut, so for simplicity we restrict our attention to this region, defining it to be 
any pixel within 30° of (260°, 0°) in Galactic coordinates. Summing all free-free emission 
in this region and dividing by the total Ha intensity in this region, we estimate that the 
ratio of radio brightness to Ha intensity kg is 9.3 ± 3.2 /iK at i^-band for the spinning 
dust fit, with similar values for the other MCMC models. The uncertainty comes from the 
variance of the ratio from pixel to pixel; increasing the signal-to-noise threshold decreases 
the uncertainty somewhat but does not significantly affect the central value. While the 
central value is consistent with the prediction of 11.4 fiK R~^ within this uncertainty, it is 
also compatible with a reduced electron temperature of 5500 K, an overestimation of the 
reddening correction by At = 0.3, or some combination of the two. 



4.4. Spinning Dust Emission 

We find that in order to best fit the 408 MHz data, the spinning dust fit from the 
five-year MCMC process needs to have its peak frequency adjusted downward by 14% from 
Usd = 4.9 GHz to Usd = 4.2 GHz, nearly independent of the offset used for the map. For 
this value, the frequency at which the flux from the spinning dust component alone peaks 
is 21 GHz. We have not found any improvement in the fit from including 'warm' spinning 
dust with a peak near 40 GHz, as found by |Dobler &: Finkbeiner 2008b The 408 MHz data 



- 16 - 



also introduces some tension, such that the spinning dust model no longer is such a large 
improvement inside the Galactic plane; in this region xl = 1-80 for the spinning dust model, 
compared to xl = 2-61 for the "base" fit, for 8.7 effective number of degrees of freedom (see 



Kunz et al. 2006 Gold et al. 2009 for detailed description of effective d.o.f.). These values 
are for the fitted offset of 7.4 K for the 408 MHz map; using a larger offset value of 11.4 
K provides slightly better fits {Axl = 0.008) for the spinning dust models, while a smaller 
offset value of 3.4K provides slightly better fits {Axl = 0.074) for models without spinning 
dust. 

The steepening synchrotron model fits the combination of WMAP and 408 MHz data 
nearly equally as well as the spinning dust model, with xl = 1-81 iii the Galactic plane. 
The amplitude of Cg is large in the Galactic plane, implying a change of spectral index 
greater than one per e-fold increase in frequency. This is a sharper change than models of 
synchrotron steepening predict from aging effects, and so the physical motivation for the 
model is unclear. 

The ARCADE data are not directly comparable to the MCMC fits, due to their greatly 
different beam and sky coverage. The spinning dust component of the fits for the two regions 
in the Galactic plane, however, does peak in flux at 22 GHz, consistent with the location 
of the MCMC peak. The relative amplitude is more difficult to ascertain. For ARCADE, 
spinning dust is 29% (at / = 33.8) or 43% (at / = 93) of the total flux at 22 GHz, but with 
the large beam it is impossible to say whether the spinning dust component is relatively 
diffuse or localized on the Galactic plane. For the MCMC fits to WMAP data, the mean 
spinning dust fraction is considerably lower, at 18% inside the KQ85y7 mask, which suggests 
that spinning dust may be patchy. Outside of the KQ85y7 mask, the MCMC fits show a 
mean level of spinning dust consistent with zero within the uncertainty of the fit. 



4.5. The Haze 



In its low frequency bands, WMAP observes an excess of emission above what was 
predicted by scaling the 408 MHz to higher frequencies using the expected spectral index 
for synchrotron emission. Determining the exact nature of this emission has proven difficult; 
WMAP has generally treated it as a hard (flatter spectrum) synchrotron component without 
attempting to explain the origin of such a component. Other suggestions have involved 



combinations of different types of spinning dust (Finkbeiner 2004 Dobler & Finkbeiner 



2008b), though there is typically still a residual "haze" even after those components are fit 



out (Dobler & Finkbeiner 2008a) 
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It has been argued that this remainder low-frequency emission has an elhpsoidal shape 
and is consistent with hard synchrotron emission, possibly from dark matter annihilation in 



the core of the Galaxy (Hooper et al. 2007). There has been tentative detection of a haze in 



gamma-rays using preliminary data from the Fermi telescope (Dobler et al. 2010). 



Interpretation of polarization information toward the center of the Galaxy is difficult, as 
depolarization through line-of-sight changes in the orientation of the magnetic field can affect 
the signal significantly. Nonetheless, we search for a hard component in the polarization data 



using a simplified version of the low- resolution MCMC fit of Dunkley et al. (2009a), shown 
in Figure |8| We do not detect any significant change of synchrotron spectral index as a 
function of Galactocentric distance. 

This special fit was done at HEALPix A^side = 16 using only WMAP polarization data, 
so as to be insensitive to any uncertainties regarding the presence or absence of spinning 
dust. The fit attempts to model the sky as a sum of three power-law foregrounds: a soft 
synchrotron component with /3 = —3.1, a hard synchrotron component with /3 = —2.39, and 
a dust component with /3 = +2.0. These power-law indices were those suggested by the work 



of Dobler fc Finkbeiner (2008a). 



The results of the fit are shown in Figure [8] Residuals after the fit are small compared 
to the noise, and over all bands the mean reduced P^r pixel is 1.1. For comparison, 
the synchrotron and dust templates used for polarization cleaning are shown in the right 
column of the figure. The MCMC result for the soft synchrotron template appears to be 
essentially a noisy version of the synchrotron template, indicating that i^-band indeed is a 
good proxy for polarized synchrotron emission. For dust, the MCMC and template results 
differ somewhat. The MCMC hard synchrotron results show no spatial structure beyond 
WMAP^s noise pattern, and are consistent with the level of noise bias expected in a map of 

P = y/Q^ + f/2. 

Figure [9] shows the frequency spectrum of polarized emission for elliptical regions around 
the Galactic center. In these regions the polarization direction is nearly vertical, and so the 
Stokes U parameter is negligible for bands K through V and small for VT-band. The spectra 
for three different regions are shown, sized 10° x 5°, 20° x 10°, and 30° x 15°. We find 
no evidence for emission other than soft synchrotron {/3 = —3.2) and dust (/3 = +2.0), in 
particular, no "haze" component appears to be necessary for polarization. 



- 18 - 



5. Foreground Systematics and Tests 

5.1. Pipeline Simulation 

A full simulation of the WMAP instrument was used to test the MCMC and template 
cleaning methods, in order to investigate the interaction of systematics in both time-domain 
data and sky maps. Starting with a set of synthetic sky inputs for the CMB and foregrounds 
(described below), the scanning of the instrument was applied to the inputs to produce a 
timestream of data, which was then put through the same entire calibration and map-making 
pipeline as used for real data. 

A random CMB realization was created, starting from the publicly available best-fit 
cosmological parameters of a ACDM model to the combination of five-year WMAP data 



with supernovae and baryon acoustic oscillations. The CAME software package (Lewis 



et al. 2000) was used to generate a model power spectrum and then synfast (Gorski et al. 



et al. 


2000 


2005 


was 1 



2005 ) was used to generate the random sky realization. 

Several foregrounds were then added, using high resolution templates. A synchrotron 



intensity template was constructed from the 408 MHz data of Haslam et al. (1981), and 
scaled to higher frequencies with a spectral index with both spatial variations and steepening, 
in order to test the effects of fitting a simpler model to complicated synchrotron spectral 
features. A free-free template was made from an extinction-corrected version of the Ha 



map of Finkbeiner (2003), with a few bright high-latitude sources removed, and assuming 



a spectral index of /3 = —2.15. The dust template is the 94 GHz prediction of model 8 of 



Finkbeiner et al. (1999), scaled to other WMAP frequencies with a spectral index of (3 = 2.0. 



Once the simulation inputs were generated, they were passed through a simulation of 
WMAP^s scan strategy, including such effects as thermal gains and baselines in the time- 
ordered data, loss imbalance and bandpass mismatches, and detector noise with a 1// com- 
ponent. This simulated time-ordered data was then processed and analyzed in exactly the 
same way as real observations. 



Figure 10 shows a comparison between the "true" simulated input sky maps and the 
output maps after the map-making process. These are used to test the template cleaning 
method, as the simulated input foregrounds are generated with structure on scales smaller 
than the templates used for cleaning. However, no effects due to residual foreground con- 



tamination are seen; the cosmological parameters used as input are recovered. Figure 11 
then compares the results of the MCMC foreground fit to the input foreground behavior. 
The largest difference found between the input and output maps from the simulation is in 
the Galactic plane. This difference is a fraction of a percent of the total intensity, and is 
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entirely consistent with the expected uncertainty in the gain reconstruction. 

The MCMC reconstructs the foregrounds to within the MCMC error, which includes 
large covariance between synchrotron and free-free brightness. The most important system- 
atic deviation was in the reconstructed synchrotron spectral shape, parameterized with Ps 
and Cs- This is largely because the simulated model spectrum was more complicated than 
a power-law with constant steepening. This resulted in a bias in the recovered Ps bias of 
approximately +0.2 in the Galactic plane. This bias was still within the MCMC errors. 



5.2. Testing Beam Systematics with Six-month Maps 

Over the course of a full year, the WMAP satellite's scan pattern is such that most 
points on the sky are observed with a nearly uniform distribution of orientations. The 
distribution is most symmetric at the ecliptic poles, and least symmetric on the ecliptic 
plane. Fortuitously, the Galactic center lies near the plane of the ecliptic, with a large angle 
between the planes of the Galaxy and the echptic. The result is that the year can be divided 
into halves, where WMAP^s scanning direction when observing the inner Galactic plane is 
rotated 180° between the two halves. 

This means that maps made from such six-month segments of data are sensitive to 
beam asymmetries, particularly those where the beam is not equal to itself rotated 180°. 



This effect is largest in i^-band. Figure 12 shows the measured difference of the beam 



between the six-month intervals, a simple beam model which recreates the effect, and sky 
maps of the residuals between six-month sky maps and a full year of observation. 

We used this effect to investigate the sensitivity of foreground fitting to beam systemat- 
ics. For the first five years of data, each year was divided into six months of one scan direction 
relative to the Galactic center, and six months where the scan direction was reversed. These 
were then stacked to produce two sets of five-year maps, where the scan directions along the 
ecliptic have the greatest relative asymmetry. The MCMC foreground fitting was then run 
for both sets of maps. 



The result is shown in Figure 13 Since the largest beam difference is in i^-band, low 
frequency foregrounds are most strongly affected. The spectral index inferred for synchrotron 
shows a small gradient across the Galactic plane, with amplitude ~ ±0.1 for \b\ < 5°. The 
effect on the CMB is limited; variance outside the KQ85y7 mask is less than 480 fxK^ (an 
order of magnitude smaller than intrinsic variance of the CMB), and most of this is from 
MCMC variations in the dust model. We also emphasize that the six-month intervals were 
chosen to maximize this asymmetry, which is not seen when full years of data are used to 
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make maps. 



6. Point Source Catalogs 



As for the five-year analysis, two separate methods have been used for identification of 
point sources from skymap data and two separate point source tables have been produced. 
The first method has been used in all WMAP data releases and is largely unchanged from 



the five-year analysis (Wright et al. 2009). The seven-year signal-to- noise ratio map in 



each wavelength band is filtered in harmonic space by hi/{bjCl + Cf°^^^), (Tegmark & 



de 01iveira-Costa||1998 Refregier et al.||2000 ) , where hi is the transfer function of the WMAP 



beam response (Jarosik et al.|2010), Ci is the CMB angular power spectrum, and Cf°^^'^ is 



the noise power. The filtering suppresses CMB and Galactic foreground fluctuations relative 
to point sources. For peaks in the filtered maps that are > 5ct in any band, the unfiltered 
temperature maps are fit with the sum of a Gaussian profile and a planar baselevel. The 
Gaussian amplitude is converted to a source flux density using the conversion factors given 



in Table 2 of Jarosik et al. (2010), and flux density uncertainty is calculated from the la 
uncertainty in the fit amplitude. The source is entered into the catalog if the fit source 
width is within a factor of two of the beam width. Flux density values are entered for bands 
where they exceed 2a. A point source catalog mask is used to exclude sources in Galactic 
plane and Magellanic cloud regions. This mask has changed from the five-year analysis. 
A map pixel is outside of the new mask if it is either outside of the diffuse component of 
the seven-year KQ85y7 temperature analysis mask or outside of the five-year point source 
catalog mask. This mask admits 82% of the sky, compared to 78% for the five-year version. 
We identify possible 5 GHz counterparts to the WMAP sources by cross-correlating with 



the GB6 (Gregory et al. 1996), PMN (Griffith et al. 1994 1995 Wright et al. 1994 1996) 



Kiihr et al. (1981), and Healey et al. (2009) catalogs. A 5 GHz source is identified as a 



counterpart if it lies within 11' of the WMAP source position (the mean WMAP source 
position uncertainty is 4', and can be twice as large for faint sources near the detection 
threshold). When two or more 5 GHz sources are within 11', the brightest is assumed to be 
the counterpart and a multiple identification flag is entered in the catalog. 

The second method of point source identification is the CMB-free method originally 



applied to one-year and three- year V- and M^-band maps by Chen & Wright (2008) and 



to five-year V- and IV-band maps by Wright et al. (2009). The method used here is that 



applied to five-year Q, V, and W maps by Chen & Wright (2009). The V- and W^-band 



maps are smoothed to Q-band resolution. A special ILC map is then formed from the three 
maps using weights such that CMB fluctuations are removed, flat-spectrum point sources are 
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retained with fluxes normalized to Q-band, and the variance of the ILC map is minimized. 
The ILC map is filtered to reduce the noise and suppress large angular scale structure. Peaks 
in the filtered map that are > 5a and outside of the seven-year point source catalog mask 
are identified as point sources, and source positions are obtained by fitting the beam profile 
plus a baseline to the filtered map for each source. Source fluxes are estimated by integrating 
the Q, V, and W temperature maps within 1°.25 of each source position, with a weighting 
function to enhance the contrast of the point source relative to background fluctuations, 
and applying a correction for Eddington bias due to noise. Detected sources were identified 
with sources in the five-year WMAP five-band catalog ( Wright et al.|[2009 ) and the five-year 



QVW catalog Chen & Wright (2009) if the positions agreed within 15'. They were also 



correlated against the 5GHz GB6, PMN, and Kiihr et al. (1981) catalogs to identify possible 



5 GHz counterparts within 15'. Optical identifications were made by searching the NASA 
Extragalactic Database. 

The seven-year five-band point source catalog is presented in Table [5] and the seven- 
year QVW point source catalog is presented in Table [6j The five-band catalog contains 471 
sources, the QVW catalog contains 417 sources, and the two catalogs have 346 sources in 
common. For comparison, the five-year five-band catalog contained 390 sources, the five-year 
QVW catalog contained 381 sources, and they had 287 sources in common. Differences in 
the source populations detected by the two search methods do not appear to be mainly due 
to spectral index differences. The distribution of spectral index in the five WMAP bands for 
the sources that are only in the five-band catalog is similar to that for the sources common 
to both catalogs. The differences are thought to be largely caused by Eddington bias in 
the five-band source detections due to CMB fluctuations and noise. At low flux levels, the 
five-band method tends to detect point sources located on positive CMB fluctuations and to 
overestimate their fluxes, and it tends to miss sources located in negative CMB fluctuations. 



This was shown by application of the method to simulated skymaps (Wright et al. 2009) 



and its effect is also seen in the comparison by Chen & Wright (2009) of five-year fluxes from 



the five-band method with those from the CMB-free method in V-, and H^-bands. 



7. Conclusions 

Even with all the uncertainty regarding foregrounds in the Galactic plane, we find no 
evidence for foreground contamination outside our current KQ85y7 analysis mask. Further, 
the cleaning process requires only three simple power-law foregrounds, and leaves no more 
than 15 /xK of residuals in the CMB temperature power spectrum. 

We find no evidence of polarized foregrounds beyond those from soft (steep-spectrum) 
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synchrotron and thermal dust emission. In particular, we see no indication of an energetic 
population of synchrotron-emitting electrons near the Galactic center. 

Additional years of data have allowed us to detect a combined 62 new point sources 
using two techniques, a 12% increase from the five-year data release. A total of 346 point 
sources are in common between the two techniques. 

More and more evidence is indicating that within a few degrees of the Galactic plane, 
the behavior of low-frequency foregrounds is complicated and has not been completely un- 
derstood. WMAP data show a rapidly steepening spectrum from 20-40 GHz, which may be 
explained as emission from spinning dust grains. The leading systematic, beam asymmetry, 
does not appear able to alter the spectrum enough to eliminate the need for spinning dust or 
a similar component. ARCADE data appear consistent with the spinning dust explanation, 
although some discrepancies remain as to the relative strength of the emission. More data 
at frequencies where spinning dust emission is expected to be strongest (10-40 GHz) would 
be very helpful. 

The WMAP mission is made possible by the support of the Science Mission Directorate 
Office at NASA Headquarters. This research was additionally supported by NASA grants 
NNG05GE76G, NNX07AL75G SOI, LTSA03-000-0090, ATPNNG04GK55G, and ADP03- 
0000-092. This research has made use of NASA's Astrophysics Data System Bibhographic 
Services. We acknowledge use of the HEALPix, GAME, and GMBFAST packages. 
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Table 1. ILC Coefficients Per Region^ 



Region i^'-band Ka-hcuid Q-band y-band W^-band 






0. 


,1495 


-0, 


.7184 


-0.3188 


2, 


.3071 


-0, 


.4195 


1 


-0. 


,0035 


-0, 


.2968 


-0.1963 


2, 


.0533 


-0, 


.5567 


2 


0. 


,0258 


-0, 


.3368 


-0.3162 


1 


.8368 


-0, 


.2096 


3 


-0. 


,0945 


0, 


.1772 


-0.6087 


1 


.5541 


-0, 


.0281 


4 


-0. 


,0771 


0, 


.0881 


-0.4149 


0, 


.9559 


0. 


,4480 


5 


0. 


,1928 


-0, 


.7451 


-0.4538 


2, 


.4673 


-0, 


.4612 


6 


-0. 


,0918 


0, 


.1946 


-0.5586 


1, 


.0227 


0. 


,4332 


7 


0. 


,1533 


-0, 


.7464 


-0.2033 


2, 


.2798 


-0, 


.4834 


8 


0. 


,2061 


-0, 


.2979 


-1.5705 


3 


.5678 


-0, 


.9056 


9 


-0. 


,0889 


-0, 


.1241 


-0.0816 


1, 


.2066 


0, 


.0880 


10 


0. 


,1701 


-0, 


.8610 


-0.1825 


2, 


.8264 


-0, 


.9530 


11 


0. 


,2358 


-0, 


.8467 


-0.6020 


2, 


.8336 


-0, 


.6206 



^The ILC temperature (in thermodynamic units) at pixel p 
of region n is T'„(p) = Yli=i Cn,iT'^ip), where C are the coeffi- 
cients above and the sum is over WMAP''s frequency bands. 
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Table 2. Template Cleaning Temperature Coefficients 



DA"" 


h 


62 (/iK R 1) 


&3 




/iff'^ (yUK R ^) 


Ql 


0.234 


1.206 


0.203 


-3.26 


7.12 


Q2 


0.232 


1.240 


0.201 


-3.30 


7.13 


VI 


0.048 


0.791 


0.466 


-3.63 


7.20 


V2 


0.045 


0.772 


0.483 


-3.64 


7.21 


Wl 


0.000 


0.436 


1.277 




7.24 


W2 


0.000 


0.430 


1.291 




7.24 


W3 


0.000 


0.438 


1.257 




7.24 


W4 


0.000 


0.432 


1.285 




7.24 



^WMAP has two differencing assemblies (DAs) for Q- 
and V^-bands and four for l^-band; the high signal-to-noise 
in total intensity allows each DA to be fitted independently. 

^Power law slope relative to i^-band, as derived from h\\ 
M^-band values are less than -4. 

*^Free-free to Ha ratio at i^-band, as derived from h\ and 
62- The expected value for an electron temperature of 8000 



K is 11.4 /iK R-^ (Bennett et al. 2003) 



Table 3. Template Cleaning Polarization Coefficients 



Band 






02'' 




Ka 


0.3202 


-3.13 


0.0144 


1.43 


Q 


0.1683 


-3.13 


0.0177 


1.54 


V 


0.0613 


-2.93 


0.0355 


1.50 


W 


0.0412 


-2.41 


0.0770 





^The coefficients are dimensionless and pro- 
duce model maps from templates. 

'^The spectral indices refer to antenna temper- 
ature. 
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Table 4. Foreground Power Spectrum Parameters 



Component Be [AiK^]"" m N [jiK] 



Synchrotron EE 


271 ± 


31 


-0.73 


± 


0.04 


0.109 


± 


0.001 


Synchrotron BB 


130 ± 


8.6 


-0.61 


± 


0.02 


0.107 


± 


0.001 


Dust EE 


17.7 ± 


2.5 


-1.13 


± 


0.06 


0.065 


± 


0.001 


Dust BB 


6.41 ± 


1.1 


-0.65 


± 


0.06 


0.066 


± 


0.001 



'^Quoted errors are only statistical uncertainty from the fit- 
ting process. 




Fig. 1. — Comparison of seven-year masks to five-year masks. At the top KQ75 and KQ75y7 
are compared, and at the bottom KQ85 and KQ85y7. Green regions are masked in both the 
seven-year and five-year masks, yellow regions are newly masked in the seven-year masks, 
and red regions are masked in the five-year masks but no longer in the seven-year masks. 
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Fig. 2. — Difference map between the seven-year ILC map and the five-year ILC map. 
Small-scale differences are consistent with pixel noise; large-scale differences are consistent 
with a change in dipole of 6.7 /iK. 
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Fig. 3. — Galactic signal component maps as determined by the Maximum Entropy Method 
(MEM) analysis. On the left are the input prior maps, and on the right are the output MEM 
maps. From top to bottom are the synchrotron, free-free, and dust components. While the 
output maps show many features of the prior at higher latitudes, there are clear differences 
in regions of strong emission. 
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Fig. 4. — Difference maps between the seven-year MEM foreground maps and the five-year 
MEM foreground maps. Apart from a small dipole shift and noise fluctuations, the only 
visible feature is a small shift of 0.17% of i^'-band flux from free-free to synchrotron. 
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(l,b) = (34°, 2°) (l,b) = (93°, 0°) 




1 10 100 1 10 100 



V (GHz) 

Fig. 5. — Galactic emission from two regions in the Galactic plane. ARCADE (triangles), 
WMAP (stars), and 408 MHz data (square) are all shown, smoothed to a common resolution. 
Upper panels show antenna temperature (absent a monopole component). The black line 
is a power-law connecting 408 MHz to 22 GHz (/3 = -2.48 for the left panel, /3 = -2.41 
for the right panel), which is divided out in the bottom panels to better show deviations 
from power-law behavior. Red lines show the result of a fit to the data using three power 
law components for foregrounds (representing synchrotron, free-free, and dust). Blue lines 
show the fit resulting when an extra component representing spinning dust is added. Solid 
lines show the total flux, with individual components shown by dashed lines (synchrotron), 
dotted lines (free-free), and dot-dashed lines (dust plus spinning dust). Errors in the data 
are dominated by systematics and highly correlated between data points, but are estimated 
to be 5 — 15%, depending on experiment. 
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Fig. 6. — Power spectrum of the difference between foreground-reduced maps. Q-band minus 
VT-band is shown here, with a point source contribution subtracted off. Note the changing 
scale between panels. Red points with error-bars are averages over bins with Ai = 50. 
Deviations from zero are below 100 fxK^ outside the KQ85y7 mask, and the upper bound to 
foreground contamination in the foreground-reduced maps is 15 /iK. 
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Fig. 7. — Power spectra of polarized foreground components as determined by the MCMC 
model. On the left are Cf^ and on the right are Cf^; for foregrounds these should be 
of comparable magnitude. The black dotted lines are the foreground fit to raw three-year 
WMAP data from Page et al. (2007), and the red dotted lines are the combined foreground 
and noise fit to MCMC maps from this work, with coefficients given in Table |4j Synchrotron 
results are in good agreement with the previous analysis. The seven-year dust results spectra 
appear to have a higher amplitude, but the signal-to-noise for £ > 10 is 2.8 or less for dust. 
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Fig. 8. — Comparison of the templates used for polarization cleaning to a low- resolution 
(iVside = 16) MCMC fit to polarization data using a three-component model with fixed spec- 
tral indexes to search for any hard synchrotron component. The left column shows the results 
of the MCMC fit to polarization data using three components: soft synchrotron (/3 = —3.1) 
at top, hard synchrotron (/3 = —2.39) at middle, and dust (/3 = +2.0) at bottom. For com- 
parison, the right column shows the templates used for polarization cleaning: synchrotron at 
top and middle, and dust at bottom. All plots are of polarization intensity P = a/Q^ + f/^, 
with a logarithmic scale from 1 to 100 /xK. Synchrotron intensity is measured at a reference 
frequency of 23 GHz, and dust intensity at 94 GHz. The MCMC maps are noisy, and have 
been corrected for a noise bias in P caused by noise in Q and U. Excess noise in the plane 
of the ecliptic due to the scan pattern is also clearly visible in the MCMC fits. Given the 
noise level, hard synchrotron emission does not appear to be significant. 
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Fig. 9. — Frequency spectrum of polarized emission around the Galactic center. Average 
antenna temperature of Stokes Q is shown for three oval regions defined by + (26)2 < 
10°, 20°, 30°, where I and b are Galactic longitude and latitude. Stokes U is negligible at 
all frequencies except VF-band. Errorbars indicate statistical uncertainty from the diagonal 
part of the pixel-pixel noise matrix. Dotted lines show the sum of a synchrotron component 
with P = —3.2 and a dust component with ^ — +2.0; in all cases this two-component model 
is sufficient to explain the observations. 
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T(mK) 

Fig. 10. — Comparison between a simulated input sky and the resulting maps after scanning 
and map-making, i^-band is shown; differences in other bands are at least four times smaller. 
The only visible structure, along the Galactic plane, is entirely consistent with residuals from 
gain reconstruction within the quoted uncertainties (0.2%). 
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Fig. 11. — Comparison between the input foreground spatial and spectral behavior and that 
recovered by the MCMC fit. Upper left: difference between MCMC result and input Tg + Tf. 
Upper right: difference between MCMC result and input T^. Lower left: difference between 
MCMC result and input (3s- Lower right: difference between MCMC result and input Cg. 
The main feature is that the simulated synchrotron model contained more steepening in the 
synchrotron spectrum than the model allowed for, which then biases the recovered Ps by 0.2 
in high signal-to-noise regions. The apparent bias off the Galactic plane only occurs where 
the signal-to-noise is low and the parameter error is larger than the bias. 
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Fig. 12. — Flight data and a simple model for differences between maps made with six 
months of data and those made with a full year. Top left: difference between observed K- 
band beam and 180° rotated i^'-band beam (scale is ±5%). Top right: difference between 
a model beam consisting of a sum of Gaussians and its 180° rotation. Middle: observed 
difference map between six months and a full year for i^-band. Bottom: simulated difference 
map created using the beam of the upper right panel. While this simple beam model does 
not completely resemble the observed beam, it qualitatively reproduces the effects observed 
in the maps. 
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Fig. 13. — Effect of beam anisotropy on the MCMC foreground fits, using stacks of six- 
montfi maps. Pixels near tfie boundary of the six-month scans are masked (gray) due to 
poor coverage. Top left: difference in MCMC synchrotron temperature. As the combination 
of synchrotron and free-free is largely constrained to match i^-band, the free-free difference 
is nearly the opposite of this map. Top right: difference in MCMC synchrotron spectral 
index. Away from the Galactic plane this map is mostly noise, but a slight gradient with 
= ±0.1 is visible near the plane. Bottom: difference in MCMC CMB temperature. 
Most of the variation is noise in the MCMC dust model, rather than due differences between 
the six-month maps. 
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